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Abstract. The ferroelectric order and magnetic field induced effects observed 
in the spiral phase of MnW04 are described theoretically. It is demonstrated 
explicitly that the Dzyaloshinskii-Moriya antisymmetric interactions contribute 
to the correlation between spins and electric dipoles in the incommensurate and 
commensurate ferroelectric phases of magnetic multiferroics. However, other 
single-site symmetric interactions are shown to be involved in the magnetoelectric 
process, suggesting the possible existence of an electric polarization originating 
from purely symmetric effects. 



1. Introduction 

A typical feature of multiferroic materials undergoing a transition to an elliptic spiral 
ferroelectric phase [110^3 , is the existence of spectacular magnetoelectric effects, such 
as the polarization flops observed in TbMnOs [4J and DyMnOs L5j , or the sign rever- 
sal of Py disclosed under magnetic field in ThMn20^ [6]. The theoretical description 
of these effects requires knowledge of the order-parameter symmetry associated with 
the ferroelectric transition. Besides, the orientation of the applied magnetic field with 
respect to the magnetic spins influences the stability range of the spiral phase and 
the polarization flop process. This property was recently illustrated by remarkable 
magnetic field induced effects observed in the ferroelectric phase of MnWO/^: The 
stability range of the phase depends on the direction of the applied field [Tl 13 [H |10] , 
which induces a high- field polarization- flop transition [11 . Here we show that these 
directional effects result from the order-parameter symmetry associated with the ferro- 
electric phase of MnW04, the stability of which depends on the respective orientations 
of the magnetic field and magnetic easy axis. 

It is of fundamental interest to understand what is the microscopic mechanism be- 
hind the magnetoelectric coupling in multiferroics. It is generally believed [3 , that 
in RMnOa (R=Tb, Dy, Gd) the antisymmetric Dzyaloshinskii-Moriya (DM) inter- 
action is the microscopic origin for the ferroelectric polarization. For RMn205 the 
spins are almost collinear in the main ferroelectric phase and it is claimed [12 that 
symmetric exchange striction induces ferroelectricity. To gain insight into the micro- 
scopic mechanism for ferroelectricity in MnW04 we express the order-parameters in 
terms of the magnetic spins and provide the correspondence between spins and elec- 
tric dipoles. This analysis confirms that the Dzyaloshinskii-Moriya (DM) interactions 
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[131 [H contribute to the electric polarization in the incommensurate spiral phase 
but, additionally, shows that they participate as well in the formation of the dipole 
moments in the commensurate ferroelectric phases of magnetic multiferroics. We also 
find that the DM interaction is not the only microscopic source of the polarization: 
Other symmetric interactions are shown to be involved in the magnetoelectric process, 
suggesting the possible existence of a polarization induced by purely symmetric effects. 



2. P^AF3^AF2^AF1 Transitions 

Below its wolframite- type paramagnetic (P) structure, of monoclinic P2/cl' symmetry, 
MnW04 undergoes three successive magnetic phase transitions [T, 'S^ at 13.5 K (Tat), 
12.7 K (T2) and 7.6 K (Ti). They lead, respectively, to an incommensurate magnetic 
phase (AF3), an incommensurate elliptical spiral phase (AF2) displaying an electric 
polarization P/ /h^ and a commensurate magnetic phase (AFl). The transition wave- 
vectors [7 are k^nc = (-0.214, ^,0.457) for AF2 and AF3, and kcom = \, \) for 
AFl. They are associated with two bidimensional irreducible corepresentations of the 
P2/cl' space group [15 , denoted by F^^ and F^^ in [7 , whose generators are given 
in Table I. 



and r*^2 ^Yve 



Table 1. Generators of the irreducible corepresentations V 
P2/cl' paramagnetic space group, deduced from the irreducible representations 
Ti and T2 of Gjz = rriy, given in Kovalevs tables [15| for \<.inc and kcom • T is the 
time reversal symmetry. For AF2 and AF3, £ ~ 0.4577r, £1 ~ 0.4187r and £2 = 2£. 
For AFl, £ = £i = 7t/2 and £2 = tt. 




The complex amplitudes transforming according to F and F , which form 
the order-parameter components, are denoted by Si = ^le^^S ^1 = Sie~^^^ and 
62 = 526*^2^ = S2e~^^'^. For k = k^^c the order-parameter invariants ^1 = S'f, 
^2 = S2 and ^3 = S1S2 cos(2(/:)), with (f = Oi — $2^ yield the Landau expansion: 

$i(T, Si,S„ ^) = $io(T) + ^5? + ^St + fsl + ^Sl + 
+y5?5|cos(2^) + ^StS^cos^{2^) + ... (1) 

Minimizing shows that five distinct stable states, denoted I-V, may arise below 
the P-phase for different equilibrium values of 81^82 and summarized in Fig. 1(a). 
The theoretical phase diagram shown in Fig. 1(b) gives the location of the phases in 
the (71,0^1 — 0^2) plane. 

Neutron diffraction data [7 show that the incommensurate magnetic AF3 phase 
induced by F^^ corresponds to = and 7^ 0, coinciding with phase I in Fig. 1(a). 
Its structure has the antiferromagnetic grey point-group 2y/myV , involving a doubling 
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Figure 1. (a) Connections between the magnetic point-groups of phases I-V 
induced by + and equihbrium conditions fulfilled by the order-parameter 
in each phase, grey rectangles indicate ferroelectric phases, (b) Phase diagram 
deduced from the minimization of $1 given in Eq.JT]), in the (7, ai — Q;2)-plane. 
The phases are separated by second-order transition curves crossing at the four- 
phase points Ni and N2. 



of the 6— lattice parameter, and an incommensurate modulation of the spin density 
in the (x,z)-plane. The AF2 spiral phase, induced by + [7J, corresponds to 
phase II in Fig. 1(a) with Sf 0^ S2 ^ 0, (p = (2n + l)f and the magnetic symmetry 

2yl'. Adding the dielectric part of the free-energy density = SPySiS2 sin(p-\- 

yy 

to Eq.([T]) yields the equilibrium polarization 

= ±5elyStSI (2) 

which changes its sign for opposite senses of the spiral configuration, as observed 
by Sagayama et al [9^. The order-parameter §2^ activated at the P^AF3 transition, 
is frozen at the AF3^AF2 transition, i.e. independent of temperature below T2. 
Therefore, Eq.(|2]) expresses a linear dependence of Py on Sf, since SSyySf acts as a 
temperature independent coupling coefficient. 

In the AF2 and AF3 phases ^1 can be truncated at the fourth degree since 
the eighth degree invariant is necessary only for stabilization of phase V (of Fig. la). 
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Putting ai = ai(T— To) and a2 = a2(T— T/v) and minimizing <l>i of Eq.([T]) with respect 
to and 5*2 one obtains the equihbrium values for the order-parameter components 
Sf and SI. 

AF3 (T2 <T<TAr)corresponds to 

nl/2 



= and Sl = ± 
For T<T2 (AF2 phase) one obtains 



(Tn - T) 



St = ±~a{T2-Ty^' and SI = ±[^{T^ - 

P2 



/2 



with 



a = { 



tt27l 



_^l^Nl/2 

v2 ) 



and T2 



a27iTN + ai/32To 



(3) 
(4) 
(5) 



/^i/^2 - 7i a27i + ai/32 

Inserting Eq.Qinto Eq.([2]) one finds that the spontaneous polarization in the 
AF2 phase varies as 

P^{T) = ±A{T2-T)^^^ (6) 

with A = S£yyd[ ^^ {T^N- T2 ) ] "'^ ^ ^ • "^^^ dielectric permittivity follows a Curie- Weiss- 
type law around T2 with 



-yy 



(T) 



AT) 



l-C 



1 - D 



with C = 



and D = C 



J-N - 1 

T-T2_ 
Tn — T2 
T2-T 



for T>T2 



for T<T2 



(7) 



Fig.2 shows that Py{T) and £yy{T) perfectly fit the experimental curves reported 
by Taniguchi et al. [8 . This confirms the hybrid character of the ferroelectricity in 
spiral magnets, recently found [T71 [18] also for TbMnOs and TbMn2 05: The square- 
root temperature dependence of the polarization (Eq.(j6])) and the Curie- Weiss type 
behaviour of £yy expressed by Eq.([7j) are typical for a proper ferroelectric transition, 
whereas the low value of Py ^ measured at 10 K [16] is of the order found 

in improper ferroelectrics. 

Let us now turn to the AFl phase. Neutron diffraction results [7 indicate that the 
transition to the commensurate AFl phase triggers a decoupling of Si and 82- Using 
Table 1 one finds that the lock-in at Ti, induced by F^^, gives rise to the additional 
invariant ^4 = cos(4^2)- The Landau expansion associated with the transition to 
AFl is: 



^2(T, 52, 02) = $2o(T) + ^Si + ^Sl + 

+ 1^4 cos(4^2) + ^Sl cos2(4^2) (8) 
4 o 

The equations of state show that three commensurate phases, denoted by F- 
to-IIF, displaying a fourfold increased unit-cell (b+c,c-b,2a+c), may appear below 
Ti. The AFl phase corresponds to phases F or IF, stable for cos(46>2) = +1 or — 1, 
respectively, both described by the magnetic space group Ca2/c. Phase IIF, stable for 

cos (4^2) = ^/f /2 ; has the symmetry CaC. 
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Figure 2. Fit (lines) of the temperature dependence of the spontaneous 
polarization Py{T) and the dielectric permittivity £yy{T) using Eq.® and 
Eq.(|3. The points are data from Taniguchi et al. [E]- The fit parameters are 
A = 21.8/iCm-2K-i/2^ c = 0.001, D = 0.0001 and = 12.3. 



3. Directional magnetic field (magnetoelectric) effects 

The AFl and AF2 order-parameters allow describing the magnetoelectric effects 
[SI El [El [m [19^ observed in MnW04. The magnetic phase diagram can be 
calculated by adding the magnetic part of the free-energy and the coupling invariants 

^1 = ^MftM - BM + i^iMfSl + i^[MfSl (9) 

to the Landau-expansion Eq.([T]) or (|8j). /i is the paramagnetic susceptibility tensor 
and i=x,y,z. 

Due to the anisotropy of the magnetic free-energy the AF2 stability range depends on 
the angle ^0 = \tan~^ \2iixz{l^zz — l^xx)~^\ between B and the magnetic easy axis in 
the paramagnetic phase. If B is at an angle ^ with the x-axis, one finds: 

T2 W - T2 = eB\l - - Vt2_4^cos2(^ - ^0)]} (10) 

where £ = {piK + 7iA^')(/3iai + 7i<^2)"\ ^ = x + l-^zz and d = ii^xl-^zz - l^lz- 
For s < the AF2 stability range is maximum for B along the easy axis = ^o)- 
It decreases when ^ increases from ^0 to ^0 + 7r/2, reducing to the stability range 
at zero field if /J^xz^fJ^xx « l-^zz- Such variation has been observed in the AF2 phase 
[TOl [20I [21], in which ^0 ^ 35^ coincides with the direction of the spins in the (x.z)- 
plane. When B is at an angle with the (x,z)-plane the AF2 stability range decreases 
when ip increases from ?/^ = 0to?/; = 7r/2, as reported experimentally JOl [20] , 
To get an overview about the magnetic field dependence of the various phases in 
MnW04 it may be useful to consider a simplified version of the free energy expansion 
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by neglecting the non-diagonal parts of jl. Then by minimizing the free energy with 
respect to Mi one obtains for the magnetization in x,y or z-direction 

and the phase transition temperatures are shifted under applied magnetic field 
Bi as 

T^(Bi) = T^iO)--^Bf (12) 

where := Ti,T2 or Tn and a > is the bare expansion coefficient of the 
second order term in Eq.([T]). Eq. (p!2]) shows, that the phase transition temperatures 
depend quadratically on the applied magnetic field and the sign of the magnetic field 
shift of Tq, depends only on the sign of the coupling coefficient tZi^ which itself can be 
determined from the changes of Mi with temperature (see Eq. (fTT]) ). Fig l3] sketches 
this behaviour. 




Figure 3. Magnetizations (a) and magnetic field dependence (b) of transition 
temperatures calculated from Eas. (|ll|) and (|12|) . 



Eqs. (pT]) and (p!2]) describe the magnetic phase diagram (see e.g. Fig. 5 of Ref [20] ) 
perfectly. The strongest downshift of Ti occurs in x-direction, which corresponds also 
to the strongest negative anomaly in M^^ at Ti, implying /^^^ >> /^^ > as seen from 
Fig. 3 of Ref [20]. This ffgure also shows that hZy < dit Ti (My displays an upwards 
anomaly at Ti) and therefore Ti{By) increases with applied magnetic ffeld. In this 
way all particular features of the phase diagram can be reproduced. 
In the following we will brieffy discuss the effect of a very high magnetic field By. Above 
a threshold field 5f , given by ^i(Ti, 5f , SI) - EfM^^ = ^2(^1, 5|) - EfM^^^', the 
AF2 phase switches to the high-field phase III', which cancels Py and gives rise to the 
polarization: 

Pi"' = -5'el,Slhin{2ei) (13) 

deduced from the dielectric free-energy $f = <5'P^5| sin(26'2) + The onset 

of P^^^ occurs in correlation with the vanishing of Py at the first-order AF2^Iir 
fiop-transition, consistent with the Py Px polarization-ffop [8 observed above 10 
T. Decreasing By below B^^ switches back phase III' to AF2. If By is canted at an 
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angle (j) with respect to the 6-axis in the (x,z)-plane, the magnetoelectric coupling 
vPxMxMySiS2 sin(/? induces the polarization 

P,^ = -elyVii-^BlSlS^^ sin ^« sin 20 (14) 

Canting oppositely By from the 6- axis {(j) ^ —(j)) reverses P^. Increasing again 
By above By^ yields an opposite sign for P^^^ in phase IIP, as observed by Taniguchi 
et al [llj . 



4. Relating the order-parameters to magnetic spins 

To gain insight into the nature of the microscopic interactions, let us express and 
^2 as a function of the magnetic spins in the commensurate phases I'-IH'. Denoting 
Si — Ss the spins associated with the eight Mn^^-ions of the corresponding fourfold 
primitive monoclinic unit-cell, one can write Si = s^a -\- s^b -\- s^c (i=l-8), where a, h 
and c are lattice vectors. Projecting on and the matrices transforming the 
5^'^'^-components gives: 



= L^'^ + iLl^' Si' = -L\ + ill 

82"^" = -Ll^' + iLl^' S2' = L\ + iL\ (15) 

where the Si^ (i=l,2 m=a,b,c) represent different forms of Si and 62. The L-^ 
(i=l-4) are projections of the vectors 



Li = Si - S2 - S7 + S8 (16) 
L2 = Si - S2 + S7 - Ss 
L3 = S3 - S4 + S5 - S6 
L4 = S3 - S4 - S5 + S6 

In phase F(S'i^'^'^ = 0,cos4^2'^'^ = 1) has s^'^ = ^2'^ = Sj'^ = 5g'^ = 0, 

h b h b r\ Ci-,c a,c a,c a,c ^ b b b b tj. ■ 

53 = S4 = 55 = = 0, S3 = —54 = —55 = Sq and Si = —S2 = —s^ = Sg. It gives 

T 1 ttu CI Ci-,b,c ^ A/)a,b,c -1 \ a,c a,c a,c a,c o,c 

In phase II (61 = 0,cos4^2 = ^1 = "-^2 = -^7 = "-^8 = ~"^3 = 
54'^ = Sg'^ = — Sg'^ and s\ = —82 = —s^ = Sg = S3 = —54 = = —Sq. It yields 

^2"'' = -Ll^" + iL^2" and ^2^ = L\ + ill 

Figs. 4(a) and 4(b) show that the magnetic structure of the AFl phase of 
MnW04 coincides with the antiferromagnetic order of phase IP (Fig. 4(b)), since 
the cancellation of the (a,c)-spin-components on the sites 1,2,7 and 8 in phase F 
(Fig. 4(a)) was not observed [7 in AFl. The lack of spin components along b reported 
in this phase [7J may be due to their relativistic origin, as suggested in [7J. In phase 
jjj,^^^a,5,c ^ Q^0a,b,c ^ o,7r/4), identified as the high-field phase of MnW04, the 

• I'-i • • a,c a,c a,c a,c a,c a,c a,c a,c 

equilibrmm spms are = —82 = 8j — —8^ , 53 = —8^ = —8^ = 8q , 
81 = —82 = 57 = — Sg and 53 = —54 = 55 = — Sg, with 62^'^ = —L^'^ + ^^2'^ and 

52' = ill 

Other commensurate structures can take place for k = kcom when the spin components 
associated with Si and S2 order simultaneously. Two among these structures, denoted 
IV' ((9f = 0,6>| = 7r/2) and V ((9f = 7r/2,6>^ = 0), display the symmetry 2y and 
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Figure 4. Calculated magnetic structures of (a) phase I' and (b) phase IP, 
described in the text. The two additional grey arrows oriented in (x,z)-plane 
located at the two atoms in the bottom left of (b) indicate the previously refined |7] 
magnetic structure of the AFl-phase. The coordinates of the eight Mn^+ ions 
corresponding to the spins si — sg, in the primitive monoclinic unit-cell of the two 



structures are for ion - 
6{i?/, I), 7(1.?/. I). 8(1.1- 



),2(i,jy,t),3(f,y,t),4(f,jy,|),5(i,jy,|) 



the same form of Py ex ^(^1 §2* — §2 §1*) = SiS2sm(f as in the AF2 phase. 
Although these structures are not stabihzed in MnW04, they are lock-in limits of 
the incommensurate AF2 structure, and can be used for investigating the microscopic 
origin of ferroelectricity in the spiral magnetic structure. Since and S2 are both 
realized by three independent combinations of spin components and taking into 
account the equilibrium relationships si = — S2, S3 = — S4, S5 = — se and S7 = — sg 
which hold for all commensurate structures of MnW04, Py is expressed as the sum of 
nine terms: 

Py = Sdsf + sf - sf - sf) + 5,{sf + sf - sf - sf) + 
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- 44) + 




- stsl + 44 - 44) + ^6(44 - 44 + 
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-44 + 44 -44) + 58(44 -44 + 




,59(44 


c a 1 c a a c\ 


(17) 



The ^1 — 84 terms in Eq.([T7j) are symmetric invariants involving a single atom. 
Their origin is entropic and due to on-site interactions. The — 5g terms repre- 
sent Dzyaloshinskii-Moriya (DM) antisymmetric coupling interactions between neigh- 
bouring pairs of spins 81,85 (or the equivalent pair 83,87), 81,83 (85,87) and 83,85 
(81,87). The DM interaction [22l [23] is currently assumed as the microscopic source 
of the polarization in the spiral structure of magnetic multiferroics [TTl [13] . Our re- 
sults confirm explicitly this view in MnW04, but show that other symmetric effects 
are also involved in the formation of the electric dipoles. Furthermore, equilibrium 
relationships of the spin components in phases IV' (^i = S7 = 0, 53'^ = 55'^ = 0) 
and V (53 = 55 = 0, 5^'^ = s^'^ = 0) preserve the symmetric and DM contributions 
in Eq. ([T7|) . This indicates that the interactions giving rise to the polarization in the 
commensurate ferroelectric phases of multiferroic compounds are of the same nature 
than in the spiral phases. The effect of the incommensurability in the x-z plane of the 
AF2 phase should result in further averaging of the spin densities without modifying 
essentially the invariants in Eq. ([T7|) . One should note that the DM interactions, as 
well as the Katsura-type contribution [13 e^j x 8^ x Sj (where e^j is the distance vector 
joining neighboring Mn atoms i and j), cancel in Eq. ([T7|) when all spin components 
are cancelled except ^i and Sj in phase IV', or 53 and 55 in phase V, Py keeping 
a finite value Py = 6i{si^ — s^'^) and Py = Si{s^'^ — sg^) in phases IV' and V, re- 
spectively. Therefore symmetry considerations predict the existence of a polarization 
induced by interactions being neither of Dzyaloshinskii-Moriya nor of Katsura-type, 
although such effects may encounter considerable restrictions at the microscopic level. 



5. Summary and conclu8ion 

In summary, the present work clarifies the nature of the ferroelectric order 
occurring in the spiral phase of MnW04, and gives a theoretical description 
of the field-induced effects observed in this compound. It confirms that the 
antisymmetric Dzyaloshinskii-Moriya interaction is involved in the formation of 
dipolar moments in the incommensurate and commensurate ferroelectric structures 
of magnetic multiferroics. It also shows that other symmetric effects participate in 
the ferroelectricity observed in these compounds, suggesting the possible existence 
of an unconventional ferroelectricity in magnetoelectric materials originating from 
purely symmetric interactions. Symmetric exchange interactions may mediate 
magnetoelectric coupling in the E-type commensurate perovskite manganese oxide 
compounds [24j. For these materials a large ferroelectric polarization was predicted 
[25, 26J and recently found in TmMn03 [24J. 

Our theoretical description of the magnetoelectric effects in MnW04 differs in 
important aspects from the theoretical approach to this compound proposed by 
Harris [27J. We determine from pure symmetry considerations the irreducible degrees 
of freedom (co-representations) involved in the observed sequence of phases, the 
corresponding order-parameter symmetries and the form of the transition free-energies 
and related magnetic, dielectric and coupling contributions, which allow describing the 
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observed magnetoelectric effects. At last we take into account the actual positions of 
the Mn ions in order to establish the connection existing between our phenomenological 
order-parameters and the magnetic spins or electric dipoles. Harris follows an opposite 
procedure, starting from the actual magnetic structure which permits construction 
of allowed spin functions. These functions are then related to the order-parameter 
components, the transition free-energy and coupling to the polarization being deduced 
from semi-empirical considerations. The advantage of our approach is that it provides 
the full set of stable states allowed by the order-parameter symmetries (one of which 
being stabilized under high magnetic field), a detailed topology of the corresponding 
phase diagram, a faithful description of the critical behaviour, including the specific 
pseudo-proper (not improper) character of the induced polarization [17 , which has 
been overlooked by Harris. It also yields a precise explanation of the various 
magnetoelectric effects observed in MnW04 (not described by Harris) and an explicit 
determination of the different types of interactions contributing to the polarization. 
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